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1 ! ABSTRACT 

Recent calculations indicate that the cohesive energy of condensed matter increases 
<— j , with magnetic field strength and becomes very significant at magnetar-like fields (e.g., 

Ck fO keV at 3 x fO 14 G for zero-pressure condensed iron). This implies that for suffi- 

q , ciently strong magnetic fields and/or low temperatures, the neutron star surface may 

?— j ■ be in a condensed state with little gas or plasma above it. Such surface condensation 

c/3 , can significantly affect the thermal emission from isolated neutron stars, and may lead 

, ^ , ■ to the formation of a charge-depleted acceleration zone ( "vacuum gap" ) in the magne- 

tosphere above the stellar polar cap. Using the latest results on the cohesive property 
04 ' of magnetic condensed matter, we quantitatively determine the conditions for surface 

condensation and vacuum gap formation in magnetic neutron stars. We find that con- 
densation can occur if the thermal energy kT of the neutron star surface is less than 
qq , about 8% of its cohesive energy Q s , and that a vacuum gap can form if f2- Bp < (i.e., 

fT) ■ the neutron star's rotation axis and magnetic moment point in opposite directions) 

qq" ! and kT is less than about 4% of Q s . For example, at B — 3 x 10 14 G, a condensed 

■ Fe surface forms when T < 10 7 K and a vacuum gap forms when T < 5 x fO 6 K. 
\ Thus, vacuum gap accelerators may exist for some neutron stars. Motivated by this 

result, we also study the physics of pair cascades in the (Ruderman-Sutherland type) 
\ vacuum gap model for photon emission by accelerating electrons and positrons due to 

V~j ■ both curvature radiation and resonant/nonresonant inverse Compton scattering. Our 

rS ' calculations of the condition of cascade-induced vacuum breakdown and the related 

■ pulsar death line/boundary generalize previous works to the superstrong field regime. 

We find that inverse Compton scatterings do not produce a sufficient number of high 
energy photons in the gap (despite the fact that resonantly upscattered photons can 
immediately produce pairs for B > 1.6 x 10 14 G) and thus do not lead to pair cascades 
for most neutron star parameters (spin and magnetic field). We discuss the implica- 
tions of our results for the recent observations of neutron star thermal radiation as 
well as for the detection/non-detection of radio emission from high-B pulsars and 
magnetars. 

Key words: radiation mechanisms: non-thermal - radiation mechanisms: thermal - 
stars: magnetic fields - stars: neutron - pulsars: general. 



1 INTRODUCTION 

Recent observations of neutron stars have provided a wealth of information on these objects, but they have also raised many 
new questions. For example, with the advent of X-ray telescopes such as Chandra and XMM-Newton, detailed observations 
of the thermal radiation from the neutron star surface have become possible. These observations show that some nearby 
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isolated neutron stars (e.g., R X J1856. 5-3754) appear to have featureless, nearly blackbody spectra ( Burwitz et al. 2003: 



an 



van Kerkwiik fe Kaplan! 120071 ) . Radiation from a bare condensed surface (where t he overlying atmosphere has negligible 



optical depth) has b een in voked to explain this nearly perfect blackbody e missi on (e.g ., Burwi tz et al J booil Mori & Rudermanl 



optical dcptli) nas been invoked to explain tins nearly pcrlcct blackbody emission (e.g., curwitz et , 
20031 : iTurolla et all I200I Ivan Adelsberg et allboost : IPerez-Azorin et aDl2006l : lllo et alJl2007t but 



Ruderman 



21)l7? f( an 



alternative v iew). How ever, whether surface condensation actually occurs depends on the cohesive properties of the surface 



matter (e.g.. lLaill2001 



Equally puzzling are the observations of anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) (see 
Woods fe Thompsonll2005l for a review) . Though these stars are believed to be magne tars, neutron stars with extremely strong 
magn etic fields (B > 10 14 G), they mostly show no pulsed radio emission (but see Camilo et al. 20061 . 2007 ; Kramer et al 



2007) and their X-ray radiation is too strong to be powered by rotational energy loss. By contra st, several high-B radio pul- 



sars with inferred surface field strengths s imilar to those of magnetars have been discovered (e.g., Kaspi fc McLaughlin 2005: 



Vranevsevic. Manchester, fe M elrose 2007). A deeper understanding of the distinction between pulsars and magnetars requires 



further investigation of the mechanisms by which pulsars and magnetars radiate and of their magnetospheres where this emis- 



matter in strong 


map 


netic fields (e.g., 


Ruderman & Sutherland 


1975; Arons & ScharlemannI 1979; Cheng & Ruderman 198C; 


Usov & Melrose! 


1996 


; Harding & Muslimovlll998l;lGil. Melikidze 


& Geppertl2003l; Muslimov & Harding 2003; Beloborodov & Thompson 


2007]). For example, 


depending on how strongly bound the surface matter is, a charge-depleted acceleration 


zone ("vacuum 



gap") above the polar cap of a pulsar may or may not form, and this will affect pulsar radio emission and other high-energy 
emission processes. 

The cohesive property of the neutron star surface matter plays a key role in these and other neutron star processes and 
observed phenomena. The cohesive energy refers to the energy required to pull an atom out of the bulk condensed matter 
at zero pressure. A related (but distinct) quantity is the electron work function, the energy required to pull out an electron. 
For magnetized neutron star surfaces the cohesive energy and work function can be many time s the corresponding terrestrial 
values, due to the strong magnetic fields thread ing the matter (e.g., Rudermanl 19741 : Lail 200 lh . 

In two recent papers |Medin fc Laill2006al lbl hereafter ML06a,b), we carried out detailed, first-principle calculations of 
the cohesive properties of H, He, C, and Fe surfaces at field strengths between B — 10 12 G to 2 x 10 15 G. The main purpose of 
this paper i s to investigate sev eral important astrophysical implications of these results (some preliminary investigations were 
reported in lMedin fc Lai 2007). This paper is organized as follows. In Section[2]we briefly summarize the key results (cohesive 
energy and work function values) of ML06a,b used in this paper. In Section [3] we examine the possible formation of a bare 
neutron star surface, which directly affects the surface thermal emission. We find that the critical temperature below which 
a phase transition to the condensed state occurs is approximately given by fcT cr it ~ 0.08Q S , where Q s is the cohesive energy 
of the surface. In Section [4] we consider the conditions for the formation of a polar vacuum gap in pulsars and magnetars. We 
find that neutron stars with rotation axis and magnetic moment given by ft ■ B p > are unable to form vacuum gaps (since 
the electrons which are required to fill the gaps can be easily supplied by the surface), but neutron stars with ft ■ B p < 
can form vacuum gaps provided that the surface temperature is less than fcT cr it ~ 0.04Q S (and that particle bombardment 
does not completely destroy the gap; see Section [6j|. In Section [5] we discuss polar gap radiation mechanisms and the pulsar 
death line/boundary in the vacuum gap model. We find that when curvature radiation is the dominant radiation mechanism 
in the gap, a pair cascade is possible for a large range of parameter space (in the P-P diagram), but when inverse Compton 
scattering (either resonant or nonresonant) is the dominant radiation mechansim, vacuum breakdown is possible for only a 
very small range of parameter values. Implications of our results for recent observations are discussed in Section [6] Some 
technical details (on our treatment of inverse Compton scattering and vacuum gap electrodynamics of oblique rotators) are 
given in two appendices. 



2 COHESIVE PROPERTIES OF CONDENSED MATTER IN STRONG MAGNETIC FIELDS 

It is well-known that the properties of matter can be drastically modified by strong magnetic fields. The natural atomic unit 
for the magnetic field strength, Bq, is set by equating the electron cyclotron energy fkj ce — h(eB/m e c) = 11.577 B12 keV, 
where B12 = £?/(10 12 G), to the characteristic atomic energy e 2 /ao = 2 x 13.6 eV (where ao is the Bohr radius): 

2 3 

B = = 2.3505 x 10 9 G. (1) 

For b — B/Bo > 1, the usual perturbative treatment of the magnetic effects on matter (e.g., Zeeman splitting of atomic 
energy levels) does not apply. Instead, the Coulomb forces act as a perturbation to the magnetic forces, and the electrons in 
an atom settle into the ground Landau level. Because of the extreme confinement of the electrons in the transverse direction 
(perpendicular to the field), the Coulomb force becomes much more effective in binding the electrons along the magnetic 
field direction. The atom attains a cylindrical structure. Moreover, it is possible for these elongated atoms to form molecular 
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chains by covalent bonding along the field direction. Interactions between the linea r chains can then lead to the formation of 



three-dimensional condensed matter IjRuderman 1974; R uder et al.lll994l ; lLai 2001 



The basic p roperties of magnetized condensed matter can be estimated using the uniform electron gas model (e.g.. 
Kadomtsevlll970h . The energy per cell of a zero-pressure condensed matter is given by 



£ 3 ~ -120 Z 9/5 B^ 5 eV, 

and the corresponding condensation density is 
p s ~ 560 A Z~ 3/5 B% 5 g cm -3 , 



(2) 



(3) 



where Z, A are the charge number and mass number of the ion fsee lLaill200ll and references therein for further refinements to 
the uniform gas model). Although this simple model gives a reasonable estimate of the binding energy for the condensed state, 
it is not adequate for determining the cohesive property of the condensed matter. The cohesive energy is the (relatively small) 
difference between the atomic ground-state energy £ a and the zero-pressure condensed matter energy £„, both increasing 
rapidly with B. Moreover, the electron Fermi energy (including both kinetic energy and Coulomb energy) in the uniform gas 
model, 



e F = (3/5Z)£ s ~ -73Z 4/5 B 2 2 /5 eV, 



(4) 



may not give a good scaling relation for the electron work function when detailed electron energy levels (bands) in the 
condensed matter are taken into account. 

There have been few quantitative st udies of infinite cha i ns and zero- pressure condensed matter in strong magnetic 
fields. Earlier variational calc ulations (e.g., Flowers_et_al. 1977|_ Mullet fl984) as well as calculations based on Thomas- Fermi 



type statistical models (e.g., Abrahams fc Shapiro 1991 ; FughDd^t~alT 1992 ). while useful in establishing scaling relations 



and providing approximate energies of the atoms and the condensed matter, are not adequate for obtaining reliable energy 
differences (cohesive energies) . Quantitative results f or the energies o f infinite chains of hydrogen molecules Hoo over a wide 



range of field strengths (5 > fi p) were presented in lLai et al.l (11992) (using the Hartree- Fock method with the plane- wave 



Relovskv fc Ruder ( 19961 ) (using density functional theory) . 



approximation; see also lLaill2001l for some results for Heoo) and in 
For heavier elements such as C and Fe, the cohesive energies of one dimensional (ID) chains h ave only been calculat ed at 
a few magnetic field stre ngths in the range of B = 10 12 -10 13 G, using Hartree-Fock models (|Neuhauser et al.lll987l ) and 
density functional theory 1 Jones! 19851) . T here were some discrep ancies between the results of these works, and some adopted 



Jones! 



si (|l986l ) 



a crude treatment for the band structure (iNeuhauser et al.lll987l) . An approximate calculation of 3D condensed matter based 
on density functional theory was presented in 



Our calculations of atoms and small molecules (ML06a) and of infinite chains and condensed matter (ML06b) are based 
on a newly developed density functional theory code. Although the Hartree-Fock method is expected to be highly accurate 
in the strong field regime, it becomes increasingly impractical for many-electron systems as the magnetic field increases, 
since more and more Landau orbitals are occupied (even though electrons remain in the ground Landau level) and keeping 
track of the direct and exchange interactions between electrons in various orbitals becomes computational rather tedious. 
Compared to previous density-functional theory calculations, we used an improved exchange-correlation function for highly 
magnetized electron gases, and we calibrated our density-functional code with previous results (when available) based on 
other methods. Most importantly, in our calculations of ID condensed matter, we treated the band structure of electrons in 
different Landau orbitals self-consistently without adopting ad-hoc simplifications. This is important for obtaining reliable 
results for the condensed matter. Since each Landau orbital has its own energy band, the number of bands that need to be 
calculated increases with Z and B, making the computation increasingly complex for superstrong magnetic field strengths 
(e.g., the number of occupied bands for Fe chains at B = 2 x 10 15 G reaches 155; see Fig. 16 of ML06b). Our density-functional 
calculations allow us to obtain the energies of atoms and small molecules and the energy of condensed matter using the same 
method, thus providing reliable cohesive energy and work function values for condensed surfaces of magnetic neutron stars. 

In ML06a, we described our calculations for various atoms and molecules in magnetic fields ranging from 10 12 G to 
2 x 10 15 G for H, He, C, and Fe, representative of the most likely neutron star surface compositions. Numerical results of the 
ground-state energies are given for Hjv (up to N = 10), He at (up to N = 8), Cjv (up to N = 5), and Fejv (up to N = 3), 
as well as for various ionized atoms. In ML06b, we described our calculations for infinite chains for H, He, C, and Fe in that 
same magnetic field range. For relatively low field strengths, chain-chain interactions play an important role in the cohesion 
of three-dimensional (3D) condensed matter. An approximate calculation of 3D condensed matter is also presented in ML06b. 
Numerical results of the ground-state and cohesive energies, as well as the electron work function and the zero-pressure 
condensed matter density, are given in ML06b for Hoo and H(3D), Heoo and He(3D), Coo and C(3D), and Feoo and Fe(3D). 

Some numerical results from ML06a,b are provided in graphical form in Figs. [1] [2] [3] and [4] (see ML06a,b for approximate 
scaling relations for different field ranges based on numerical fits). Figure Q] shows the cohesive energies of condensed matter, 
Q s = £\ — £ s , and the molecular energy differences, A£n ~ £n/N — £%, for He, Fig. [2]for C, and Fig. [3]for Fe; here £\ is the 
atomic ground-state energy, £n is the ground-state energy of the Hejv, Cat, or Fejv molecule, and £ s is the energy per cell 
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of the zero-pressure 3D condensed matter. Some relevant ionization energies for the atoms are also shown. Figure [4] shows 
the electron work functions (f> for condensed He, C, and Fe as a function of the field strength. We see that the work function 
increases much more slowly with B compared to the simple free electron gas model [see Eq. @] , and the dependence on Z is 
also weak. The results summarized here will be used in Section and Section U below. 



3 CONDENSATION OF NEUTRON STAR SURFACES IN STRONG MAGNETIC FIELDS 

As seen from Figs. [T] [5] and[3] the cohesive energies of condensed matter increase with magnetic field. We therefore expect that 
for sufficiently strong magnetic fields, there exists a critical temperature Tcxix below which a first-order phase transition occurs 



between the co ndensate and the gaseous vapor. This has been investigated in detail for hydrogen surfaces fsee lLai fe Salpeter 



1997; Lai 2001), but not for other surface compositions. Here we consider the possibilies of such phase transitions of He, C, 
and Fe surfaces. 

A precise calculation of the critical temperature T cr i t is difficult. We can determine T cr it approximately by considering the 
equilibrium between the condensed phase (labeled "s") and the gaseous phase (labeled "g") in the ultrahigh field regime (where 
phase separation exists). The gaseous phase consists of a mixture of free electrons and bound ions, atoms, and molecules. 
Phase equilibrium requires the temperature, pressure and the chemical potentials of different species to satisfy the conditions 
(here we consider Fe as an example; He and C are similar) 

P s = P g = [2n(Fe + ) + 3n(Fe 2+ ) + ■ • ■ + n(Fe) + n(Fe 2 ) + n(Fe 3 ) + ■ • -]kT , (5) 
Us = fie + M(Fe + ) = 2ne + M(Fe 2+ ) = • • • = ^(Fe) = ^(Fe 2 ) = ^(Fe 3 ) = ■ • ■ , (6) 
where we treat the gaseous phase as an ideal gas. The chemical potential of the condensed phase is given by 
fi a — £ a + P S V S ^ £ s ,o , (7) 
where £ s is the energy per cell of the condensate and £ a ,o is the energy per cell at zero-pressure (we will label this simply as 
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£ s ). We have assumed that the vapor pressure is sufficiently small so that the deviation from the zero-pressure state of the 
condensate is small; this is justified when the saturation vapor pressure P sat is much less than the critical pressure P cr i t for 
phase separation, or when the temperature is less than the critical temperature by a factor of a few. 
For nondegenerate electrons in a strong magnetic field the number density is related to fj, e by 

1 Mc /fcT f-n L fiw ce \ f°° dp z ( -p\ \ . . 

e"«/ > , g nL cxp — — / — exp — — (8) 



2npf ^»n L ^y kT J J h ^\2m e kT 

n L =0 

1 ^/ fcT tanh-' lT ' x - 



l—e^ kT , (10) 



2TYpf ) \ T , 



2np 2 \T, 

where g nL = 1 for i%l = and g„ L = 2 for ul > are the Landau degeneracies, Ar e = (2nh 2 /m e kT) 1 ' 12 is the electron 
thermal wavelength, and the last equality applies for kT <JC hcu ce . The magnetic field length is po = (hc/eB) 1 ^ 2 . For atomic, 
ionic, or molecular Fe the number density is given by 



1 



^ f d*K exo 



^exp(-^^)z int (Fc A ), (12) 



TA 

with the internal partition function 
A£ A , 



Z int (Fe A ) = ^exp(-^i) . (13) 

i 

and A£a,% = £a,% —£a- Here, the subscript A represents the atomic, ionic, or molecular species whose number density we are 
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100 p-TTT! 




10 100 1000 

B (10 12 G) 

Figure 3. Cohesive energy Q s = £i — E s and molecular energy difference Afjv = £n /N — £i for iron as a function of the magnetic field 
strength. The symbol Qoo represents the cohesive energy of a one-dimensional chain, and I\ and I2 are the first and second ionization 
energies of the Fe atom. Below 5 X 10 12 G, our results for and Q 3 become unreliable as Qoo and Q 3 become very small and approach 
numerical errors for £jy and £ s . 



calculating (e.g., Fe2 or Fe + ) and the sum J^. is over all excited states of that species. Also, \re = (2-kH 2 /MAkT) 1 ^ 2 is the 
Fe particle's thermal wavelength, where Ma = NAM is the total mass of the particle (N is the number of "atoms" in the 
molecule, A is the atomic mass number, and M — m p + m e ). The vector K represents the center-of-mass momentum of the 
particle. Note that we have assumed here that the FeA particle moves across the field freely; this is a good approximation for 
large Ma- The internal partition function Zi nt represents the effect of all excited states of the species on the total density; in 
this work we will use the approximation that this factor is the same for all species, and we will estimate the magnitude of 
this factor later in this section. 

The equilibrium condition fi s = /u(Fe) for the process Fe Sj00 + Fe = Fe Sj00+ i yields the atomic density in the saturated 
vapor: 

^H^H exp {-kf) z ^' (14) 

where Q B = £1 — £ s is the cohesive energy of the condensed Fe. The condition Nfj, s = ^(Fejv) for the process Fe Sj0 o + Fejv = 
Fe s ,tx>+iv yields the molecular density in the vapor: 

,„ , / NAMkT\ 3 / 2 ( S N \ 
where 

S n =Sn- N£ s = N[Q S - (£1 - £ N /N)] (16) 

is the "surface energy" and £n/N is the energy per ion in the molecule. The equilibrium condition /i(Fe n+ ) = /i e +/i(Fe' n+1 - )+ ) 
for the process e + Fe n+ = Fe' n+1 - )+ , where Fe n+ is the nth ionized state of Fe, yields the vapor densities for the ions: 
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Figure 4. Numerical result for the electron work function as a function of the magnetic field strength, for He, C, and Fe infinite chains. 



n(Fe ; 



2 +i 



2-kcLq 



lm e kT ( 

V^ exp v 



kTJ 



n(Fe 4 



(18) 



and so on. Here, b = B/Bq and ao is the Bohr radius, and 7 n = f( n _i)+ — £ n + represents the ionization energy of the nth 
ionized state of Fe (i.e., the amount of energy required to remove the nth electron from the atom when the first n — 1 electrons 
have already been removed). The total electron density in the saturated vapor is 



n(Fe + ) + 2n(Fe 2 



+ • 



(19) 



The number densities of electrons [Eq. (|19[) ] and ions [e.g., Eqs. (|17[) and (118[l ] must be found self-consistently, for all ion 
species that contribute significantly to the total vapor density. The total mass density in the vapor is calculated from the 
number densities of all of the species discussed above, using the formula 

p 3 = AM [n(Fe) + 2n(Fe 2 ) + ■ ■ ■ + n(Fe+) + n(Fe 2+ ) + ■ ■ •] . (20) 

Figure [5] (for Fe) and Fig. [6] (for C) show the the densities of different atomic/molecular species in the saturated vapor 
in phase equilibirum with the condensed matter for different temperatures and field strengths. These are computed using 
the values of £n/N, £„, and £ n + presented in ML06a,b and depicted in Figs. [2] and [3] As expected, for sufficiently low 
temperatures, the total gas density in the vapor is much smaller than the condensation density, and thus phase separation is 
achieved. The critical temperature T cr it, below which phase separation between the condensate and the gaseous vapor occurs, 
is determined by the condition p s = p g . We find that for Fe: 



3 x 10 6 , 10 7 , 2 x 10 7 K for Bi 



5, 10, 100, 500, 1000, 



for C: 

T crit ~ 9 x 10 4 , 3 x 10 5 , 3 x 10 6 , 2 x 10 7 K for B 12 = 1, 10, 100, 1000 
and for He: 

T crit ~ 8 x 10 4 , 3 x 10 5 , 2 x 10 6 , 9 x 10 6 K for B 12 = 1, 10, 100, 1000 



(21) 
(22) 
(23) 
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In terms of the cohesive energy, these results can be approximated by 

kT clit ~0.08Q s . (24) 

Note that in our calculations for the iron vapor density at B12 = 5-500 we have estimated the magnitude of the internal 
partition function factor Zint; the modified total density curves are marked on these figures as u p g x Zi nt " . To estimate Zi n t 
we use Eq. ((13} with a cutoff to the summation above some energy. For B12 = 5, 10, 100, and 500 we calculate or interpolate 
the energies for all excited states of atomic Fe with energy below this cutoff, in order to find Zi n t. The energy cutoff is 
necessary because the highly excited states become unbound (ionized) due to finite pressure and should not be included in 
Z^t (otherwise Zint would diverge) . In principle, the cutoff is determined by requiring the effective size of the excited state 
to be smaller than the inter-particle space in the gas, which in turn depends on density. In practice, we choose the cutoff 
such that the highest excited state has a binding energy \£a,i\ significantly smaller than the ground-state binding energy \£a\ 
(typically 30% of it). As an approximation, we also assume that the internal partitions for Fejv molecules and ions have the 
same Z m t as the Fe atom. Despite the crudeness of our calculation of Zj nt , we see from Fig. [S] that the resulting T cr i t is only 
reduced by a few tens of a percent from the T cr it value assuming Zint — 1 ■ 

We note that our calculation of the saturated vapor density is very uncertain around T ~ T cr i t , since Eqs. (|I4[I - (|f 8[) are 
derived for p g <C p 3 while the critical temperature of the saturated vapor density is found by setting p s = p g . However, since 
the vapor density decreases rapidly as T decreases, when the temperature is below T cr it/2 (for example), the vapor density 
becomes much less than the condensation density and phase transition is unavoidable. When the temperature drops below a 
fraction of T cr it, the vapor density becomes so low that the optical depth of the vapor is negligible and the outermost layer 
of the neutron star then consists of condensed matter. The radia tive properties of such con densed phase surfaces have been 
studied using a simplified treatment of the condensed matter (see Ivan Adelsberg et al.ll2005l and references therein). 



4 POLAR VACUUM GAP ACCELERATORS IN PULSARS AND MAGNETARS 



A rotating, magnetized neutron star is surrounded by a magnetosphere filled with plasma. The plasma is assumed to be an 
excellent conductor, such that the charged pa rticles move to screen out any electric field parallel to the local magnetic field. 
The corresponding charge density is given by |Goldreich fe Julianllf 969T ) 



pa j 



tt ■ B 

2nc 



(25) 



where ft is the rotation rate of the neutron star. 

The Goldreich- Julian density assumes that charged particles are always available. This may not be satisfied everywhere 
in the magnetosphere. In particular, charged particles traveling outward along the open field lines originating from the polar 
cap region of the neutron star will escape beyond the light cylinder. To maintain the required magnetosphere charge density 
these particles have to be replenished by the stellar surface. If the surface temperature and cohesive strength are such that the 
required particles are tightly bound to the stellar surface, those regions of the polar cap through w hich the charged particles 
are e s caping will not be replenished. A vacuum gap wi l l then develop just above the polar cap (e.g., 



Ruderman fe Sutherland 



arc escaping will not be replenished. A vacuum gap will men develop just above tne polar cap (e.g., Kudcrman & butncrlant 
i97l Icheng fe Rudermanllf98ol : lUsov fe Melroselll996l : Izhang. Harding, fe Muslimovll200ol : iGil. Melikidze. fe GeppertH2003T ) 



In this vacuum gap zone the parallel electric field is no longer screened and particles are accelerated across the gap until vacuum 
breakdown (via pair cascade) shorts out the gap. Such an acceleration region can have an important effect on neutron star 
emission processes. We note that in the absence of a vacuum gap, a polar gap acceleration zone based on space-c harge-limited 
flow may still develop (e.g., Arons fe Scharlemann 19791 ; Harding fe Muslimov 19981 : Muslimov fe Harding! 20031) . 

In this section we determine the conditions required for the vacuum gap to exist using our results summarized in Section^ 
The cohesive energy and electron work function of the condensed neutron star surface are obviously the key factors. We examine 
the physics of particle emission from condensed surface in more detail than considered previously. 



4.1 Particle Emission From Condensed Neutron Star Surfaces 

We assume that the NS surface is in the condensed state, i.e., the surface temperature T is less than the critical temperature 
T cr it for phase separation (see Section [3]). (If T > T cr i t , the surface will be in gaseous phase and a vacuum gap will not form.) 
We shall see that in order for the surface not to emit too large a flux of charges to the magnetosphere (a necessary condition 
for the vacuum gap to exist), an even lower surface temperature will be required. 



4-1.1 Electron Emission 

For neutron stars with fi-Bp > 0, where B p is the magnetic field at the polar cap, the Goldreich- Julian ch arge density is nega- 
tive at the polar cap, thus surface electron emission (often called thermionic emission in solid state physics; I Ashcroft fe Mermin 
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Figure 5. The mass densities of various atomic/ionic/molecular species and the total density (p g ) of the vapor in phase equilibrium 
with the condensed iron surface. The five panels are for different field strengths, B12 = 5, 10, 100, 500, 1000. The horizontal lines give 
the densities of the condensed phase, p s . All the vapor density curves are calculated assuming Z; nt = 1, except for the curve marked 
by "p g X Zi nt " , for which the total vapor density is calculated taking into account the nontrivial internal partition functions of various 
species. The critical temperature T cr i t for phase separation is set by the condition p g = p s . 



1976) is relevant. Let T e be the number flux of electrons emitted from the neutron star surface. The emitted electrons are 
accelerated to relativistic speed quickly, and thus the steady-state charge density is p e — — eJF e /c. For the vacuum gap to 
exist, we require \p e \ < pg.i|- (If \eJ- e /c\ > |pgj|> the charges will be rearranged so that the charge density equals Pgj-) 

To calculate the electron emission flux from the condensed surface, we assume that these electrons behave like a free 
electron gas in a metal, where the energy barrier they must overcome is the work function of the metal. In a strong magnetic 
field, the electron flux is given by 

f (26) 

where p m i n = y/2m e \Uo\, Uo is the potential energy of the electrons in the metal, e = p1/(2m e ) is the electron kinetic energy, 
and 
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Figure 6. The mass densities of various atomic/ionic/molecular species and the total density (p g ) of the vapor in phase equilibrium 
with the condensed carbon surface. The four panels are for different field strengths, B12 = 1, 10, 100, 1000. The horizontal lines give the 
densities of the condensed phase, p a . All the vapor density curves are calculated assuming Z[ n ^ — X. The critical temperature 7crit for 
phase separation is set by the condition p g = p s . 



M = e («-»i)/*T + 1 ( 27 ) 

is the Fermi-Dirac distribution function with p' e the electron chemical potential (excluding potential energy) . Integrating this 
expression gives 

where <$>= \Uo\~ p' e is the work function of the condensed matter and the second equality assumes <j> kT. The steady-state 
charge density supplied by the surface is then 

p e = --Fe = Pgj exp (C e - <j>/kT) , (29) 
c 

with 

(e kT \ 
~ o — 1 — 21 f ) — 31 + In (PoTe) ~ 30 , (30) 

where T 6 = T/(10 6 K) and Pq is the spin period in units of 1 s. For a typical set of pulsar parameters (e.g., Po = 1 and 
Te = 0.5) C e ~ 30, but C e can range from 23 for millisecond pulsars to 35 for some magnetars. Note that the requirement 
(f> ^> kT is automatically satined here when |p e | is less than |pgj|- The electron work function was calculated in ML06b and 
is depicted in Fig. [4] 



4-1.2 Ion Emission 



For neutron stars with ft ■ B p < 0, the Goldreich-Juliam charge above the polar cap is positive, so we are interested in ion 
emission from the surface. Unlike the electrons, which form a relatively free-moving gas within the condensed matter, the ions 
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are bound to their lattice sites0 To escape from the surface, the ions must satisfy three conditions. First, they must be located 
on the surface of the lattice. Ions below the surface will encounter too much resistance in trying to move through another 
ion's cell. Second, they must have enough energy to escape as unbound ions. This binding energy that must be overcome will 
be labeled £b- Third, they must be thermally activated. The energy in the lattice is mostly transferred by conduction, so the 
ions must wait until they are bumped by atoms below to gain enough energy to escape. 

Consider the emission of ions with charge Z n e from the neutron star surface (e.g., Fe + would have Z n = 1). The rate 
of collisions between any two ions in the lattice is approximately equal to the lattice vibration frequency Vi, which can be 
estimated from 

where Q p — (A-kZ 2 e^Uilm^^ is the ion plasma (angular) frequency and uj C i — ZeB / '(rriic) is the ion cyclotron frequency 
(rrii — Arrip). Not all collisions will lead to ejection of ions from the surface, since an energy barrier Sb must be overcome. 
Thus each surface ion has an effective emission rate of order 

X = v,e~ £BlkT . (32) 

The energy barrier Sb for ejecting ions of charge Z n e is equivalent to the energy required to release a ne utral atom from the 
surface and ionize it, minus the energy gained by returning the electron to the surface (e.g., Tsong 199p| ). Thus 



£b = Qs + Y^ J ' - Z ™<t> . ( 33 ) 

where Q s > is the cohesive energy, Ii > is the ith ionization energy of the atom (so that ^ — i ^ * s ^ ne ener gy required 
to remove Z n electrons from the atom), and (j> > is the electron work function. The surface density of ions is riiVi, where Ti 
is the mean spacing between ions in the solid. Thus the emission flux of Z n -ions is 

Ti = v % n % r % e- £BlkT . (34) 

The steady-state Z n -ion number density supplied by the surface is then 
Z„e 

Pi = —^—Ti = pgj exp(Ci - £ B /kT) , (35) 
c 

with 

Znevimn 



d = In 



CpGJ 



~ M + \n ^Z n ZA~ 1/2 n i 2 / s 2 (r l /a )Bi 2 1 Poy / l + 5.2 x 10~ s A- 1 Bf 2 n^ ~ 27-33, (36) 

where n2s = n;/(10 28 cm -3 ). For a typical set of pulsar parameters (e.g., B12 = 1 and Pq = 1) C, ~ 27, but d can be as 
large as 33 for magnetars with B\2 = 1000 and Pq = 8. 

All the quantities in £b were calculated in ML06b (see Figs. [2] and [3}. We find that the emission of singly-ionized atoms 
(Z n = 1) is most efficient, as £b is signficantly lower for Z n = 1 than for Z n > 1 (y) . " Ii grows much faster with Z n than 
Z n (j) does). 



4-1.3 Effect of Electric Field on Charge Emission 

The discussion in Sections 14.1.11 and 14.1.21 includes only thermal emission of charged particles from the condensed surface. A 
strong electric field, of order E s ~ £lBR/c, may be present. Since this electric field is much less than the characteristic field 
~ e/r 2 inside the condensed matter (where n is the mean particle separation), this field cannot directly rip charges off the 
surface. Nevertheless, the electric field may enhance the thermal emission of charge particles. We now estimate the magnitude 
of this effect. 

In the presence of a vacuum gap, the electric field E s at the stellar surface points outward (E s > 0) for stars with 
tt ■ Bp < and inward (E s < 0) for stars with tt ■ B p > 0. A charge Q moved to some small height z above the surface gains 
a potential energy given by U — —Q 2 /(Az) — QE s z, where the first term is due to the interaction between the charge and the 
perfectly conducting metal surface, and the second term is due to the external field0 The potential reaches a maximum value 

1 The freezing condition is easily satisfied for condensed matter of heavy elements (see Ivan Adelsberg et al.ll2005l) . 

2 In the vacuum gap, the electric field is not exactly uniform, but since the maximum U is attained at a rather small height compared 
to the gap thickness, this nonuniformity is unimportant for our consideration here. 
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t/ ma x = -|Q| 3/2 |£ s | 1/2 (37) 

at the height z — |Q/4_E S | 1 ^ 2 . Thus, compared to the E s = case, the energy barrier for particle emission is now reduced by 
the amount U max . 

Combining this consideration with the results of Sections 14.1,11 and 14,1.21 we find that steady-state charge density due to 
electron surface emission (for fl ■ B p > stars) is (cf. Jessner et*lrt] 200lf ) 



p e = PG j exp[C e -(<f>- e 3/2 \E s \ 1/2 )/kT], (38) 
and the steady-state charge density due to ion surface emission (for fl ■ B p < stars) is 

= pajexp[d - (£b - (Z n e) 3/2 \E s \ 1/2 )/kT]. (39) 
For E a ~ flBR/c, we have e 3 ^ 2 \E a \ 1 ^ 2 ~ 10 eV. This is typically much smaller than either <j> or £g. 

4.2 Conditions for Gap Formation 

No vacuum gap will form if the electrons or ions are able to fill the magnetosphere region above the polar cap with the 
required Goldreich- Julian density; i.e., the vacuum gap will cease to exist when p e — pcj or p, ; = pcj- From Eqs. (|39|l and 
(|38|) we can see that no polar gap will form if 

<f>-e 3/2 \E s \ 1/2 < C e kT~3T 6 keV (40) 
for a negative polar magnetosphere (fl ■ B p > 0), and 

£ B - (Z n ef /2 \E S \ 1/2 < dkT ~ 3T 6 keV (41) 

for a positive polar magnetosphere (fl ■ B p < 0). [For the exact expressions for C e and d see Eqs. (I30|l and 1)360 .] 

For neutron stars in general, the electron work function <j> is much less than C e kT ~ 3Te keV (see Fig. U}, so electrons 
can easily escape from the condensed surface. N o gap forms for a nega tive polar magneto sphere under neutron star surface 



conditions. (This is contrary to the conclusions of Usov fe Melrose! 19961 and Gil et al. 20031 .) The ion binding energy £b [given 



by Eq. (|33[) ]. on the other hand, can be larger than dkT ~ 3T^ keV under certain neutron star surface conditions (see Figs.[T] 
[2] and Ions can tightly bind to the condensed surface and a polar gap can form under these conditions. Figure [7] shows 
the critical temperature (determined by £b = CikT) below which a vacuum gap can form for the Fe, C, and He surfaces. 



VACUUM GAP ACCLERATORS: PAIR CASCADES AND THE PULSAR DEATH 
LINE /BOUNDARY 



Pair cascading in the magnetosphere of a pulsar is an essential ingredient for its radio emission (e.g.. iMelrosd 120041 ) . The 
pair cascade involves: (a) acceleration of primary particles by an electric field parallel to the magnetic field; (b) gamma ray 
emission by the accelerated particles moving along the magnetic field lines (either by curvature radiation or inverse Compton 
upscattering of surface photons); (c) photon decay into pairs as the angle between the photon and the field line becomes suffi- 
ciently large. To initiate the cascade an acceleration region is required; the charact eristics of this particle accelerator determine 
whether pulsar emission can operate or not (t he so-called "pulsar death line"; e.g., Ruderman fe Sutherland] 1975 : Aronsllioool : 



Zhang et alJboOCl : iHibschman fe Aronsl 120011 ). Depending on the boundary condition at the neutron star surface, there are 
two types of polar gap accelerators: If charged particles are strongly bound to the neutron star s urface by cohesive forces, 



a vac uum gap develops directly above the surface, with height h much less than the stellar radius (jRuderman fe Sutherland 
19751 ): if charged particles can be freely extracte d from the surface, a more e xtended space-charge-limited-flow (SCLF) type 



accelerator develops due t o field line curvatures ( Arons fe Scharlemann 19791 ) and the relativistic frame dragging effect (e.g. 
Muslimov fe Tsvganlll992l). B ecause the cohesive strength of matter at B ~ 10 G was thought to be negligible (ba sed on 



the result of 



Ncuhauscr et al 



200d : iMuslimov fe Harding||2003l . 12004) 



1987), m ost theoretical works in recent years have focused on the SCLF models (e.g., Aronsl 



Our results in Section [4] show that for sufficiently strong magnetic fields and/or low surface temperatures, a vacuum 
gap accelerator can form. Such a vacuum gap may be part icularly relevant for the so- c alled high-B radio pulsars, which have 
inferred magnetic fields similar to those of magnetars (e.g., Kaspi fe McLaughlin 2005 : Burgav et al. 2006). In this section we 
discuss the conditions under which a vacuum gap will be an effective generator of pulsar emission. As discussed in Section [4] 
since electrons are weakly bound to the condensed stellar surface, such a vacuum gap is possible only for pulsars with fJB p < 
(as suggested in the original Ruderman-Sutherland model). 

Our analysis is similar to the original Ruderman-Sutherland model, except that we extend our discussion of the cascade 
phys i cs to the magnetar fie ld regime, which introduces some corrections to previous works (e.g., Ruderman fe Sutherland] 
19751 : lUsov fe Melrose! fl996l ). We also consider photon emission due to inverse C ompton scattering, in addition to curvature 
radiation, in the cascade (cf. Zhang et al. 1997 . 2000l : Hibschman fe Arons 2001 ). 
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Figure 7. The condition for the formation of a vacuum gap above condensed helium, carbon, and iron neutron star surfaces, when the 
magnetosphere is positive over the poles (fi ■ B p < 0). 



5.1 Acceleration Potential 

When the temperature drops below the critical value given in Section 01 the charge density above the polar cap decreases 
quickly below pgj, and a vacuum gap results. In the vacuum region just above the surface (0 ^ z <C R), the parallel electric 
field satisfies the equation dE\\/dz ~ —4:npcj- The height of the gap h (<C R) is determined by vacuum breakdown due to 
pair cascade, which shorts out the electric field above the gap (i.e., E\\ = for z ^ h). Thus the electric field in the gap is 

E^^(h-z), (42) 

where B p = bdB p is the actual magnetic field at the pole, and differs from the dipole field Bp by a factor bd ^ 1. The potential 
drop across the gap is then 

A $ = 9*Z h ' = bd ^l h 2 . (43) 
c c 

With this potential drop, the electrons and positrons can be accelerated to a gamma factor 

7m = ; = 5.43 x 10%/13-Po = 1-23 x \^bdB^h\P^ , (44) 

where f3 Q = B p /Bq (with B Q = m 2 e c 3 /eh = 4.414 x 10 13 G the QED field), B 12 = B^/(W 12 G), h 3 = h/{10 3 cm) and P is 
the spin period in units of 1 s. The voltage drop across the gap can be no larger than the voltage drop across the polar cap 
region A$ maa; ~ (QB p /2c)(r p+ ) 2 = (QB p /2c)(r p+ ) 2 , where r p+ = r% + /b d ' 2 is the radius of the polar cap through which a 
net postive current flows: 



(§f»(25f. 

Thus the gap height is limited from above by 



d 



r d 



/W ~ -Z±= = 7.54 x 10 3 b~ 1/2 P~ 1/2 cm, (46) 



I2bd 
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where we have adopted R — 10 km. 

The above equations are for an aligned rotator. For an oblique rotator (where the magnetic dipole axis is inclined relative 
to the rotation axis), the voltage drop across the polar cap region is larger, of order (£lB p /2c)Rr p+ . But as discussed in 
Appendix A, the acceleration potential across the vacuum gap is still limited from above by A& max ~ (QB p /2c)r 2 . 



5.2 Requirements for Gap Breakdown 

There are two requirements for the breakdown of a vacuum gap. First, the photons must be able to create electron-positron 
pairs within the gap, i.e., the mean free path of photon pair-production is less than the gap height: 

l ph < h . (47) 

Second, the electrons and positrons must be accelerated over the gap potential and produce at least several photons within the 
gap. If on average only one photon is emitted with the required energy for each electron-positron pair, for instance, then the 
number of charged particles produced in the gap will grow very slowly and the gap will not break down completely. Therefore, 
we must have 

iVph > A , (48) 
where N p ^ is the number of photons emitted within the gap by each electron or positron, and A is a number of order 1-10. 



5.3 Pair Production 



The threshold of pair production for a photon with energy e is 
e 



■ sin 9 > 1 . 



(49) 



2m e c 2 

where 9 is the angle of intersection of the photon and the magnetic field. Suppose a photon is emitted at an angle 9 e . After 
the photon travels a distance z, the intersection angle will grow as z/lZ c , where 1Z C is the local radius of curvature of the 
polar magnetic field line. Thus the typical intersection angle (for a photon crossing the entire gap) is 

„ „ h 



+ i 



(50) 

10 8 P 1//2 cm, but a more complex field topology at the 



For a pure dipole field, the curvature radius is of order (Rc/Q,) 1 ^ 2 
polar cap could reduce TZ C to as small as the stellar radius. 

In the weak-field regime, when the thresho ld condition is well-satisfied (so that the pairs are produced in highly excited 
Landau levels), the mean free path is given by ( Erber 19661 ) 



4.4a / 4 

p Q sin 8 6XP I 3x 



with 



X = 



2m e c 2 



/3q sin( 



(51) 



wher e ao = h 2 / (m e c 2 ) is the Bohr radius. The condition i p h < h implies \ ^ 1/15 for typical parameters ( Ruderman fc Sutherland! 
19751 ). For stronger magnetic fields (/3q > 0.1 — 0.2), the pairs tend to be pr oduced at lower Landau levels. Using the gen- 
eral expression for the pair production rate (e.g.. iDauehertv fc Hardin dfl983h . one can check that if the threshold condition 
Eq. (149 \ is satisfied, the pair-production optical depth across the gap would also be greater than unity [for (3q — 0.1, the 



optical depth r is unity when e/(2m e c 2 ) sin 9 > 1.05, and by (3q = 0.2, r = 1 when e/(2m e c 2 ) sin 9 > 1 + 10 7 .] Thus for 
arbitrary field strengths, the condition l p h < h leads to the constraint: 

1 



2m e c 2 



(— 



> ^ 1 



15/Sq). 



(52) 



5.4 Photon Emission Multiplicity and the Pulsar Death Line/Boundary 

There several possible photon emission mechanisms operating in the vacuum breakdown, each leading to a different "death 
line" , or more precisely, "death boundary" . We consider them separately. 



5.4-1 Curvature Radiation (CR) 

The characteristic energy of a photon emitted through curvature radiation is e ~ (3/2)-f 3 hc/TZ c = 4.74 x 10 9 /3q/i3P _3 7?. ( ^ 1 eV, 
where TZq — 7?. c /(10 6 cm), and we have used 7 ~ j m [Eq. ((44}]. The emission angle is 9 e ~ 7 -1 , which is typically much less 
than h/TZ c (this can be easily checked a posteriori). Equation ((52} then reduces to 
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h > /w ph = 546P 3 ^r 



1/7 

cm. (53) 



]V D|1 ,^,^J = 17.6/^Po-X" 1 • (54) 



The rate of energy loss of an electron or positron emitting curvature radiation is Pcr = 2e 2 7 4 /(3c 3 )(c 2 /7?. c ) 2 , thus the 
number of photons emitted through curvature radiation by a single electron or positron across the gap is 

Vph ~ ~~T~c~ ~ 9hcTl c 
The condition jVp h > A [Eq. (gSJ] then gives 

h > ftxoin.e = 384A 1/3 /3~ 1/3 P 1/3 ^g /3 cm. (55) 

Thus the minimum gap height required for vacuum breakdown is h ~ max(/i m i n!P h, /i m in,e)- Combining Eqs. 1)46 p . ()53)l . and 
((551) . we have 

maxf/imin^ph, /i m in,e) < Vax ■ (56) 

This gives a necessary condition for pulsar emission and defines the pulsar "death line". For all relevant parameter regimes, 
ftmio.ph > ftmin,e, and Eq. (156)1 simply becomes /i m j niP h < /i max - The critical pulsar spin period is then 

P crit = 1.64&y 13 B? 2 /13 7^ 4/13 (l + 15/3 Q )- 2/13 s, (57) 

w here the dipole polar field i s B 12 = 2.0(P As) 1/2 , with As = P/(10~ 15 s s" 1 ). For f3 Q < 1/15 this is the same as the result 
of lRuderman fc Sutherland 1 19751 ). 

In Fig. (8j we show the death lines determined from Eq. ()56)l for the cases of TZa = 1 and TZq — WOP^ 2 (pure dipole field 
at the polar cap), with bd — 1. 

5.4-2 Resonant Inverse Compton Scattering (RICS) 

Here the high-energy photons in the cascade are produced by Compton up scatterings of thermal photons from the neutron 
star surface. Resonant scattering in strong mag netic fields (e.g. jHeroldlll979h can be thought of as resonant absorption (where 



the electron makes a transition from the ground Landau level to the first excited level) followed by radiative decay. Resonance 
occurs when the photon energy in the electron rest frame satisfies e ~ e c = h(eB /m e c) = /3gm e c 2 . The resonant photon 
energy (in the "lab" frame) before scattering is ei = e c /[y(l — cos ft)], where ft is the incident angle (the angle between 
the incident photon momentum and the electron velocity). After absorbing a photon, the electron Lorentz factor drops to 
7 e = 7/ (1 + 2/3q) 1,/2 , and th en radiatively decays isotropically in its rest frame. The characteristic photon energy after resonant 



scattering is therefore fe.g.. iBeloborodov fc Thompson 2007) 




6 = 7(1 1 I m e c 2 , (58) 

with typical emission angle ft ~ l/7e- The condition Z p h < h [see Eq. 1)520 ] becomes 

1 \ . ( h ,/TTW Q \ . „ 1 m 

For /3q > 4 this condition is automatically satisfied, i.e., resonant ICS photons pair produce almost immediately upon being 
upscattered. For /3q < 4, Eq. (|59|) puts a constraint on the gap height h. As we shall see below, most of the scatterings in the 
gap are done by electrons/positrons with 7 ~ min(7 c ,7 m ), where 7 C = e c /kT (with T the surface blackbody temperature) 
and 7 m is the Lorentz factor of a fully-accelerated electron or positron [Eq. (|44[)]. For 7 = 7 m , Eq. ([59} yields 

h > ^L.ph = 56.9P 1/3 7^ /3 /(/3 Q ) 1 '' 3 cm, (60) 

where 

For 7 = 7c we have 

h Z ^!n 2 in, P h = 169ft 6 T 6 /(/3 Q ) cm. (62) 
Combining Eqs. ()60p and ()62[1 . we find that the condition i p h < h leads to 

h > Vm, P h = max(ftW pll , ft^L.ph) • ( 63 ) 
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Figure 8. Pulsar death lines/boundaries for the CR and resonant ICS gap breakdown mechanisms. For curvature radiation, the lower 
line is for a magnetic field radius of curvature comparable to the stellar radius (1Z C ~ R) and the upper line is for a radius curvature 

1/2 

given by the dipole formula (TZq = IOOPq ). For RICS, the large "box" is for A = 1 and the small box is for A = 2; both boxes are for 
a surface temperature of 5 X 10 6 K. The unspecified neutron star parameters are taken to be unity (i.e., we set = 1 and for RICS 
1Z.§ = !)• The CR mechanism operates (and the pulsar is alive) above and to the left of the lines, and the RICS mechanism operates 
within the boxes. Radio/X-ray pulsars (ATNF catalog, http://www.atnf.csiro.au/research/pulsar/psrcatl are labeled by crosses, while 
magnetars (McGill catalog, http://www.physics.mcgill.ca/~pulsar/magnetar/main.html) are labeled by solid circles and the two radio 
magnetars are labeled by solid triangles. 



The resonant cross section for inverse Compton scattering, in the rest frame of the electron before scattering, is 
i e 2 h 



2tt z — - 5(e' - e c ) 



(64) 



where e' ~ 7^. This cross section is a ppropriate even for B p > Bq, since the resonant condition e' = e c holds regardless of 
field strength (cf. Gonthier et al. 2000l ) . The ambient spectral photon number number density near the polar cap is 



dn ph _ e*/(Tr 2 h 3 c 3 ) 



dc, 



s ti/kT _ I 



(65) 



For concreteness, consider a positron produced at z = with initial Lorentz factor 7 = 1 and accelerated to 7 = 7 m after 
crossing the full gap@ Neglecting the radiation reaction (see later), we have 7- 1 = 2(7™ - l)(z/h - z 2 /2h 2 ). The number of 



3 We can also consider the general situation where a positron (electron) is created at some location within the gap with initial Lorentz 
factor much less than -y m , travels upwards (downwards) across the gap and get accelerated to a final Lorentz factor of order 7 m . This 
would give similar result for N p ^. 
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photons upscattered through resonant ICS by the positron is given by (see Appendix B) 

f-h />oo j 

iVph ~ I dz I dei . ph a. 



o 



dc, 

d 7 I 7 — 1 



(7m - 1) ao J 1 7 3 (e E <=/ fcT ^ - 1) 

j_ / fcr \ 2 h r c 

j m \m e c 2 ) a J x 



-1/2 



( e * -1)(1- aj m / x ) V2 

where we have used 7 m 3> 1 and 
£ c 7c , nn w 1n -3 



(66) 



— = 1.09 x 10" 3 /i3 2 P T f r 1 . (67) 



Imkl 7 

Note that the second equality of Eq. (|66|l gives 



^^wV^-ra-^-r''*. (68) 

From this equation we see that for -y c — e c /kT < 7™,, dA r p h/dln7 peaks at 7 ~ 7 e , with (dAr p h/dln7) 7=7c ~ N p h, while for 
'Jc < ■y < 7™,, dN ph /d\n-y is of order (7 c /7)iVph; for 7 C > 7 m , dN ph /dlwy ~ ( r r/^ m )N ph peaks at 7 ~ 7 m . Therefore, most 
of the scatterings in the gap are done by electrons/positrons with 7 ~ min(7 c ,7 m ). Since we are interested in the regime 
e c /kT 3> 1, the integral in the last equality of Eq. (|66|) depends only on x m , and for our purpose it can be approximated as 
(n 2 /6)x m (e Xm — 1)~ . This approximation reproduces the exact integral in the x m —* limit. Thus we have 

3/2 

iV ph ,re B 4.89 x 10' 2 Pq 1 T^ f 2 P^ 2 F(x m ) , with F{x m ) = J™_ 1 . (69) 

The function F(x m ) peaks at x m = 0.874 with F maj: = 0.585. Thus the condition iVph > A necessarily requires 2.86 x 
10- 2 /3q 1 T 6 5/2 P 1/2 > A, or 

0g < /3 Q ,crit = 2.86 x 10- 2 A- 1 T 6 5/2 P 1/2 . (70) 

For a given /3q < /?Q, C rit, the condition N ph > A is equivalent to F(x) > 0.588/3q / Pq >C iU,, which limits x m to the range 
x a < x m < Xb, where x a ^ are determining by solving F(x m ) = 0.588/3Q//3Q, cr it. This condition then translates to the 
constraint on h: 

^min.c *C h *C /lraax,c , (^^-) 

where 

Vin.c = 33i- 1/2 P 1/2 T»r 1/2 cm , Vax,c = 33^ 1/2 P 1/2 T 6 - 1/2 cm. (72) 

In summary, vacuum breakdown involving RICS requires 
Pq < Aj.crit and max(/i min , p h, /imin.c) < min(/i max , /i max , c ), (73) 



where /?Q, C rit, La, ft m m, P h, ftmin.c, /imax.o are given by Eqs. J70J, (]46j), (J63J) (note that /i m in, P h = for Pq > 4), and (J75J, 
respectively. In Fig. [8] we show the pulsar death boundary when RICS is most important for initiating a cascade in the vaccum 
gap, for the cases A = 1 and A = 2, with ba — 1, 1Z& = 1, and T$ = 5. Note that in Fig. [8] we have not plotted RICS death 
boundaries for the case of a dipole radius of curvature = IOOPq^ 2 ) or a surface temperature Ta < 1; there are no regions 
of the P-P diagram where vacuum gap pair cascades are possible under these conditions. 

The pulsar death boundary depicted in Fig. [8] can be understood as follows: (i) a) The condition h^ in ph < /i max gives 

(la) P < 352 b- 3/5 1l- 2/5 f(f3 Q )- 2/5 s, (74) 

where /(/3q) is given by Eq. I|6ip. This is shown as the long-dashed line labeled (la) in Fig. [9] b) The condition h^ m ph < /i max 
gives 

(lb) P < 1.99 x 10 3 b^TZe 2 T h - 2 f(l3 Q )- 2 s. (75) 

This is shown as the short-dashed line labeled (lb) in Fig. [9] This set of conditions, (la) and (lb), is the usual requirement 
that photons emitted by an accelerated electron or positron in the gap must be able initiate pair production, (ii) a) For 
Pq < /3<2,crit, we have x a ~ 0.342 {Pq / Pq, c m) 2 , and the condition ft£j nph < ftmax.e then yields 

(Ha) P > 210 X 3/2 Tll /2 Te 3 p 3 Q /2 f(P Q ) 1/2 s. (76) 
This is shown as the dotted line labeled (Ha) in Fig. [9] b) The condition h^ in ph < /i max , c yields 
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(lib) P> 105 XTleT^fiPQ) s. (77) 

This is shown as the dot-long-dashed line labeled (lib) in Fig. [9] This set of conditions, (Ha) and (lib), together with 
Pq ?S /3Q,crit, come from the requirement for efficient photon emission by RICS in the gap. (iii) The condition /i m i n ,o < Lai 
gives 

1/2^1/2^1/2 o ^ _ n a7A ,„ fee 

Pq 

This condition is shown as the dot-short-dashed line labeled (III) in Fig. [5] (iv) The condition (5q > /3Q, C rit gives Eq. (I70|l and 
is shown as the light solid line labeled (IV) in Fig. 

Previous studies of the the pulsar death conditions for vacuum gaps where RICS is the dominant photon emission 
mec hanism have found that the RICS mechanism can lead to gap breakdown for a wide range of neutron star parameters (see, 
Zhang et alB oOO). This is contrary to our results, which show (see Figs. [8] and [9} that RICS is not a good mechanism 



(III) P < 2286" ' T;' xl> s, with x b ~ 0.874 + In ^ . (78) 



e.g 



for gap breakdown, except under very specific conditions (e.g., high surface temperatures and long rotation periods). The 
discrepancy arises because previous works did not calculate/estimate N p h (the number of high energy photons produced as a 
positron/electron crosses the gap) correctly. For example, it was implicitly assumed that photon product ion continues acros s 
the entire gap at the same rate as it does when 7 ~ y c (i.e., at the point of maximum RICS power loss) ( Zhang et al. 2000l ) . 



This assumption is invalid for 7 > y c , as is discussed above: dN p h/d\ny grows with increasing gamma factor until 7 ~ 7 C , 
and then it decreases [see Eq. (|68[) ]: therefore, dN p h/dy (which is directly related to the photon production rate -/V p h) drops 
faster than 7 -1 above 7 ~ 7 C . 

Note that the accelerating positron/electron is not radiation-reaction limited at 7 ~ y c , since the power loss due to RICS 
is significantly smaller than the power gain due to traversal across the potential drop. The power loss due to RICS is given by 

3 dn 

<ki —j^- a' res (e - a) (79) 

~ 2 ^ c d 1 ^ (SO) 

- a [ ^T+2fa ) 7 2 {e^ kT ~< - 1) " 1 ' 

At the point of maximum RICS power loss (when 7 = 7 C = e c /kT) 

P loss ( 7 = 7c ) ~ ^(l-—J == \(llL)\e-l)- 1 me c 2 (81) 
ao \ y/l + 2/3 Q / Kmc 2 J 



fcf. lDermerl (|l990l )). The power gain due to acceleration across the gap is given by 



(82) 



eS» c = ^^-(h-z)m e c 2 . (83) 



aao 

Thus 



~ 170 ( PoVPq ( 1 - 1 ] . (84) 
VlOOcnJ PQ y ^l + 2fa) 

For most pulsar parameters, P ga m ^> -Pioss [e.g., in order for 7 to reach 7 C the gap height must be at least h = 33P 1 ^ 2 T 6 ~ 1//2 cm; 
see Eq. (|67|l with x m = 1]. Therefore, there is no reason why 7 should remain near y c , the point of maximum RICS photon 
emission, as was assumed in some earlier papers. 



5.4-3 Nonresonant Inverse Compton Scattering (NRICS) 

The characteristic energy of a photon Compton-upscattered by an electron or positron of Lorenz factor 7 is e ~ ye' /(l + x), 
where x — e'/mc 2 , e ~ yet, and et is the initial seed photon energy; the pitch angle of the scattered photon is of order 
Be ~ (1 + x)/y. In the vacuum gap, most the scatterings are by electrons/positrons with 7 ~ y m on seed photons with initial 
energy a ~ 2.82fcT (see below). Substituting 

7>" £ m -.v. e m 2.82kTym /or\ 

with x m = = (85) 



l+x m mc- 
into Eq. (|52|l (which results from the requirement Z p h < h), we find 

ym ( \ a ( ^ 1 1 I «£m \ ^ ft , 1 

~2 



\l+x m J \TZ C y m I 15 
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Figure 9. The pulsar death boundaries when the resonant ICS mechanism is most important for initiating a cascade, plotted as a 
function of the dimensionless magnetic field strength /3q = B/Bq and the period P. The boundaries are shown for surface temperature 
T = 5 X 10 6 K and parameter A = 1 (the largest, bold, enclosed region), T = 5 X 10 6 K and A = 2 (the mid-sized enclosed region) and 
T = 2 X 10 6 K and A = 1 (the smallest enclosed region). The critical lines defining the edges of the region for Tg = 5, A = 1 are also 
shown. Each critical line (I)— (IV) is determined by one of Eqs. 1741 — 1781 and 1701 1. as discussed in Section 15.4.21 



Using Eq. 
0.0415/3, 



this becomes 



-1/2 D l/2^,-3/2^-l x 



5/2 



+ X m > 2 1 



(87) 



l + Xr 

The gap height is related to x m by 

h = 19.7a;^ /2 /3Q 1/2 P 1/2 r 6 " 1/2 cm. (88) 
The solution to Eq. (|87|l yields x m > a; m i n , and thus the constraint on the gap height from / p h < h is 

h > Vin, P h = 19.7zV 2 n /3g 1/2 P 1/2 T 6 - 1/2 cm. (89) 
The nonresonant part of the ICS cross section, in the rest frame of the electron before scattering, is approximately given 



by 

cr'(e') ~ ctt 



/kn(x) = ot 



/kn(i) ■ 



e +e c J \ x + Pq 

where ax is the Thomson cross-section, x = e' /(mc 2 ), and 



/kn(i) 



l + x I 2a:(l+x) \ 1 1 + 3z 

^T2^ - ln(1 + 2x) \ + Tx ln(1 + 2x) (TT2^ 



(90) 



(91) 
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is the Klein-Nishina suppression factor [/kn — 1 — 2a; for i C 1 and /kn — (S/8x)(\n 2x + 1/2) for i>l]. This agrees well 
with the calculated NR cross sections in strong magnetic fields (e.g., Gonthier et al. 2000h . 
The number of scatterings per unit length by an electron or positron is 



3 



dz _^ *^-'(7^)~0.24(-j a'(2.82kTj), (92) 

where in the second equality we have used the fact that dn p h/dmei peaks at = 2.82kT, while cr'(e') varies more slowly 
with e'. Similar to Section [5.4.21 consider a positron produced at z — with initial Lorentz factor 7 = 1 and accelerated to 
7 = 7 m after crossing the full gap. The number of scatterings produced by the positron is given by 

27-7! v/1 - 7/7™ dz 

Clearly, most of the scatterings are by positrons/electrons with 7 ~ j m , producing photons with energy 2.82fcT7^/(l + x m ) 
[see Eq. (|85|l ]. The number of photons scattered by 7 = (0.7-l)7 m electrons/positrons is 

h fdN uh \ /kT\ 3 



P h 



where 



1.3 x 10- 4 /3- 1/2 T 6 5/2 P 1/2 F(x m ,/3 Q ), (94) 



5/2 



F{x m ,Po) = 7 , „ x 2 /knM- (95) 

(x m + /^q) 2 

Now consider the vacuum breakdown condition N v h > A. For a given {3q, the function F(x m , Pq) has a maximum 
fmax(/3Q) (this maximum is approximately achieved at x m ~ 2.24 + 3/3q). Then > A requires 

P > Porit(/3Q) = 5.7 x 10 7 A 2 T 6 - 5 /3 Q F max (/3,3)- 2 s. (96) 

When this is satisfied, we additionally require 
1 -1/2 



F{x m ,f3 Q ) 

-Fmax(/9Q) 



P 



Per* (Pq) 



(97) 



which yields the solution x a < x m < X),. In terms of the gap height, we have 

^min.c *^ h> *C ^max,c , (98) 

where 

ft mi n,c = 19.7a;I /2 ^- 1/2 P 1/2 r- 1/2 cm, Vax.o = 19.7xy 2 £- 1/2 P 1/2 T- 1/2 cm. (99) 

When the neutron star surface temperature Tr, ^ 5 there are no values of /3q or P for which NRICS can initiate a cascade in 
the vacuum gap. (Only when > 9 are there any P values which permit an NRICS-initiated cascade, and even at these 
high temperatures the allowed range of (3q and P values is very small and atypical of neutron stars.) Therefore, no pulsar 
death boundaries appear for the NRICS process in Fig. [8] 



6 DISCUSSION 



It is well known that a strong magnetic field increases the binding energy of individual atom and that of the zero-pressure 
condensed matter. Very approximately, for B Bo [see Eq. (1)], the former increases as (InB) 2 while the latter scales as 
B 0A . Therefore one expects that the outermost layer of a neutron star may be in the condensed state when the magnetic field 
B is sufficiently strong and/or the surface temperature T is sufficiently low. Exactly under what conditions this occurs is an 
important question that entails quantitative calcula tions. In this paper, using our recent results on the cohesive properties of 
magnetized condensed matter I Medin fc La i 2006a b), we have established quantitatively the parameter regime (in B and T 
space) for which surface condensation occurs. Our calculations showed that there are a range of neutron star magnetic field 
strengths and surface temperatures where the condensed surface will have an important effect on radiation from these stars. For 
example, if the surface composition is Fe, then strong-field neutron stars (B > 10 13 G) with moderate (T < 10 6 K) surface 
temperatures should have atmospheres/ vapors that are effectively transparent to thermal radiation, so that the emission 
becomes that from a bare condensed surface. This may explain the nearly blackbody-like radiation spectrum observe d from 
the nearby isolated neutron star RX J1856. 5-3754 (e.g., Burwitz et al. 20031 ; van Adelsberg et al. 20051 ; Ho et al. 2007 ). 

We have also examined the conditions for the formation of a vacuum acceleration gap above the polar cap region of 
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the neutron star. The inner acceleration gap model, first developed bv lRuderman fc Sutherland! 1 19751 ). has provided a useful 
framework to understand numerous observations of radio pulsars . Most notably, the model naturally explains the phe nomenon 
of drifting subpulses observed in many pulsars (e.g., Backer|[l976l ; Deshpande fc Rankinfl999l ; Weltevrede et al. 2006h in terms 
of th e E x B circulation of plasma filaments pr o duced by vacuum discharges. Partially screened gaps have also been studied 
(e.g., ICheng fc Rudermarjfl98(il ; iGil et al-lEooil . 2006). However, it has long been recognized that the original Ruderman & 
Sutherland model is problematic since the dipole magnetic field inferred from P, P may not be strong enough to inhibit charge 
emission from the surface. Our calculations described in this paper quantify the condition for vacuum gap formation (see 
Fig.0. While this condition (i.e., T is smaller than a critical value which depends on B and composition) may not be satisfied 
for most pulsars (unless one invokes surface magnetic fields much stronger than that inferred from P, P; see Gil et al. 2006 
and references therein), it could well be satisfied for some neutron stars. In particular, the recently d iscovered high-B radio 
pulsars, having dipole surf ace magnetic fields in excess of 10 14 G and temperature about 10 6 K fe.g.. lKaspi fc Gavriill 2004; 
Kaspi fc McLaughlinll2005h . may operate a vacuum gap accelerator. On the other hand, while magnetars have similar magnetic 



field strengths, their surface temperatures are about five times larger than those of high-B radio pulsars, and therefore may 
not have a vacuum gap. In this regard, it is interesting to note that most magnetars do not show radio emission (though this 
may be because the radio pulse is beamed away from us or the because their magnetosphere plasma "overwhelms" the radio 
pulses), and the two recently detected radio magnetars have rather different radio emission properties (e.g., the spectrum 
extends to high frequency and the radiation shows high degrees of linear polarization) compared to "normal" radio pulsars. 
We may therefore speculate that a key difference between magnetars and high-B radio pulsars is their difference in surface 
temperature. In any case, our gap formation condition (Fig.O suggests that the radio emission property of neutron stars may 
depend not only on the magnetic field and rotation rate, but also on the surface temperature. 

We note that our calculation of the requirements for vacuum gap formation assumes idealized conditions. A real n eutron 
star polar cap may be immersed in a strong radiation field and suffer bombardment from high energy particles (e.g., Aronsl 



198ll : iBeloborodov fc Thompsonll2007l ). The effective cohesive energy of the surface may be somewhat smaller than what we 



used in our paper due to surface defects (Arons 2007, private communication). Whether t he vacuum gap s urviv es in realistic 
situations is unclear. It has been suggested that a partially screened gap is formed instead jGil et al.ll2003l . bood ). With small 
modifications [e.g., the potential drop given by Eq. (|43[1 is reduced], our discussion of pair cascades in the vacuum gap can be 
easily generalized to the case of a partially screened gap. 

A major part of our paper is devoted to the pair cascade physics in the vacuum gap (Section[5}. We find that pair cascade 
initiated by curvature radiation can account for most pulsars in the P-P diagram, but significant field line curvature near the 
stellar surface is needed. Although such field curvature is possible for some pulsars, it is unlikely to occur for all of them. For a 
pure dipole magnetic fiel d, only about half of all pulsars can be explained by a curvature radiation-initiated cascade. Contrary 



to previous works (e.g., Zhang et al. 2000), we find that inverse Compton scatterings (resonant or not) are not efficient in 



producing vacuum breakdown via pair cascade. 

The recent detection of the radio emission from two AXPs ijCamilo et al.ll2006l . 2007) is of great interest. The emission 
appears to be triggered by X-ray outbursts of usually quiescent magnetars. This may be due to a rearrangement of the surface 
magnetic field, which made pair cascades possible. We note that the occurrence of pair cascades depends strongly on the 
field line geometry/curvature; our study of pair cascades in the context of vacuum gap accelerators (Section [5| serves as an 
illustration of this point. 
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APPENDIX A: MAXIMUM POTENTIAL DROP FOR AN OBLIQUE ROTATOR 

For an oblique rotator, with the magnetic inclination angle a much larger than the polar cap angular size r p /R, the voltage 
drop across the polar cap is of order (£lB p /c)Rr p sin a, which is a factor of R/r p larger than the aligned case. Here we show 
explictly that the maximum potential drop across the height of the vaccum gap is still of order (tlB p /c)r p . 

We will be working in the "lab" frame, where the star is rotating. For simplicity we approximate the vacuum gap to be a 
cylinder of radius r p and height h <C R; see Fig. IA1I In reality the bases of the cylinder are not exactly circular for an oblique 
rotator, but this does not affect our conclusion. The gap is small compared to the stellar radius and we can treat it locally, 
using a Cartesian coordinate system: z along the gap height and x,y for the distance from the pole (with x being principally 
along 9 and y along <j>; 6 points in the direction from the rotational pole to the magnetic pole). The magnetic field in the 
cylinder is approximately uniform, B = —B p z. 

The electric potential inside the cylindrical gap satisfies V 2 $ = 0. The potential at the base and on the walls of the 
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Figure Al. A schematic diagram showing the polar gap structure in the cylindrical approximation. 



cylinder can be found from E = —c~[(fl x r) x B], At the top of the cylinder, the potential satisfies E z = —d&/dz — 0. 
With these boundary conditions the potential within the cylinder is completely determined. 

Without solving the complete potential problem, here we only consider the potential at the top of the cylinder in the 
limit of h 3> r p (but h <C R); this corresponds to the maximum potential drop across the gap. In this limit, the electric field 
on the top satisfies not only E z = but also dE z /dz = 0. Thus the Laplace equation below the top of the cylinder becomes 



dx 2 dy 2 



To lowest order in x/R and y/R the electric field below the base of the cylinder is given by 



E 



(n 



x B 



sin a + cos a 



R 



V » 

x + cos a— y 
R ' . 



(Al) 



(A2) 



The potential at the base of the cylinder is therefore (using E = 
the pole) 



x I x 2 iy 2 \ 

sin a— + cos a H I 

R \2R 2 2R 2 J 



-V$ and renormalizing such that the potential is zero at 



(A3) 



Since E z — on the cylindrical wall, the potential on the wall is also given by Eq. <)A3[) . The potential at the top of the 
cylinder must solve Eq. (|Aip and match the potential on the wall along the upper edge. For a circular polar cap boundary, 



given by x 2 + y 2 



the potential at the top is then 



$top = 



Q.B V R 



sin a— + cos a 
R 



2R 2 



(A4) 



From Eqs. (IA3|) and (|A4[) . we find that at the magnetic pole, |3>to P — "J'bascl — (OB p rp/2c) cos a, which is the value of 



the aligned case multiplied by cos a. Figure IA2] compares the potential at the base and top of the cylindrical gap along the x 
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Figure A2. Potential along the x (~ 8) direction, through the magnetic pole, both at the stellar surface and at the top of the vacuum 
gap, for an oblique rotator. The magnetic inclination angle is chosen to be a = n/8. The potential is measured in units of £lB p r 2 /2c, 
the value of the maximum potential drop for an aligned rotator case. 



axis. We see that although there is a large potential drop across the polar cap, the potential difference between the top and 
the base is smaller. 

Alternatively, we may examine the problem in the rotating frame, in which the potential inside the vacuum gap satisfies 
the equation V 2< 3? = 4-7rpGJ, where pa is the Goldreich- Julian charge density. The potential at the base and on the wall of 
the cylinder is $ = (since the electric field is zero there). At the top of the cylinder, we have dQ/dz — 0. These boundary 
conditions completely determine the potential inside the cylinder. For /i < r p , we expect the potential drop along the z-axis, 
A$, to grow as h 2 . But when h becomes larger than r p , the potential drop A$ will saturate to {Q,B v r 2 /2c) cos a, similar to 
the aligned case. 



APPENDIX B: SCATTERING RATE CALCULATION 

In the neutron star rest frame ("lab" frame), the electron (positron) is embedded in a radiation field with specific intensity 
I H (yii). In the electron rest frame, the radiation intensity is 

l[,((l')=(L] / ei (4), (Bl) 



where e' and £j are related by a Lorentz transformation: e' = £i7(l — /3cos0j). For a photon coming in along the Cl' direction, 
the total scattering cross section is a' = f dfK ( -fe- ] , which in general depends on Cl' and e'. The scattering rate in 

the electron rest frame is 

w = hf^- (B2) 

In the lab frame the scattering rate is dN/dt = r y~ 1 {dN /dt') (e.g., Rvbicki fc Lightman 1979t ). Using dft' /dfti = (e^/e') 2 and 
Eq. ()B1|I we have 



^ = Jd^j dei (1 - (3 cos 6i)^-a'. (B3) 

Neglecting the angle dependence of a' and assuming that the radiation field I ti is isotropic, this becomes 
dN f dn ph , 
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which is the same as Eq. (|66[1 . 
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